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SUMMARY 

This  report  considers  matters  reiated  to  the  application  of  instrument  compatibility 
checking  techniques  to  flight  test  data.  A  previously  developed  Maximum  Likelihood 
program  has  been  used  to  study  the  effects  of  the  presence  of  scale  errors,  accelerometer 
offsets  and  measurement  time  lags  using  simulated  data.  Some  additional  information 
on  the  effects  of  noise  levels  has  also  been  obtained.  The  results  have  led  to  a  suggested 
method  for  determination  of  centre  of  gravity  hcation  from  flight  data.  The  effects  of 
measurement  lags  have  been  .shown  to  have  a  motor  influence  on  extracted  instrument 
parameters  and  a  systematic  procedure  for  the  determination  df  relative  phases  has  been 
devised  and  applied  successfully  to  .simulated  data.  These  techniques  have  also  been  applied 
to  flight  data  from  a  roller-coaster  manoeuvre  and  a  set  of  relative  lag  values  clearly 
identified.  The  question  of  the  accuracies  of  extracted  instrument  parameters  and  their 
dependence  on  the  relative  lags  will  Iw  treated  more  fuliv  in  a  .subsequent  publication. 
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Sensitivity  matrix,  equation  (11). 

Linear  accelerations  in  x,  z  directions,  m/s*. 

Offset  biases  in  ax,  az,  q  etc.  measurements. 
Gravitational  acceleration,  m/s*. 

Altitude,  m. 

Time  index. 

Cost  functional  to  be  minimized,  equation  (8). 
Measurement  noise  vector,  equation  (7). 

Random  noise  in  ax,  az,  q  measurements. 

Number  of  time  points. 

Pitch  rate,  rad/s. 

Measurement  noise  matrix, 
time,  sec. 

Input  vector. 

Velocities  in  x,  z  directions. 

Airspeed,  m/s. 

State  vector. 

Estimated  state  vector. 

Reference  body  axes. 

Body  axis  coordinates  of  angje  of  attack  probe,  m. 
Acoeleiomcter  positions  rdative  to  centre  of  gravity. 
Observation  vector. 

Calculated  observation  vector. 

Angle  of  attack,  radians. 

Increment. 

Scale  factor  error  in  ax,  az,  q  etc.  measurements. 
Pitch  attitude,  radians. 

Parameter  vector. 

Standard  deviation. 

Gradient  with  respect  to  parameter  vector,  (. 
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1.  INTRODUCTION 


At  ARL,  TKe  Aircraft  Behaviour  Studies-Fixed  Wing  Group  has  for  some  years  been 
active  in  the  application  of  system  identification  methodology  to  the  extraction  of  aerodynamic 
information  from  flight  test  data.  The  group  has  acquired  or  developed  a  number  of  parameter 
and  state  estimation  programs  (Ref.  I)  and  acquired  experience  in  their  use  and  application 
to  flight  test  data  analysis. 

The  problem  of  extracting  aerodynamic  information  from  dynamic  flight  test  data,  especially 
when  the  form  of  the  aerodynamic  model  is  uncotain,  is  compounded  by  errors  in  the  measure¬ 
ments,  particularly  systematic  bias  and  scale  factor  errors.  One  approach  taken  seeks  to  remove 
instrumentation  errors  and  inconsistencies  prior  to  further  analysis,  by  making  use  of  redundant 
information  available  in  the  quantities  measured.  This  is  done  by  suitably  formulating  the  exact 
aircraft  kinematic  equations,  relating  accelerations,  velocities  and  displacements  and  including 
instrument  biases  and  scale  factors  as  unknown  parameters.  Solutions  to  the  resulting  non-linear 
state  and  parameter  estimation  problem  have  been  obtained  with  simulated  data  using  both 
Extended  Kalman  Filter  (Ref.  2)  and  Maximum  Likelihood  (Ref.  3)  methods  and  a  comparison 
of  the  two  approaches  is  documented  in  Reference4.  Both  techniques  were  shown  to  give  acceptable 
results  with  simulated  data  provided  certain  requirements  on  manoeuvre  shape,  noise  levels, 
data  length  and  sampling  rate  are  met.  Further  development  of  the  Extended  Kalman  Filter 
Approach  and  related  techniques  has  been  pursued  under  a  Research  Agreement  with  the  Uni¬ 
versity  of  Newcastle  and  reported  in  References  5  and  6.  The  present  report  concerns  some  further 
results  obtained  using  the  Maximum  Likelihood  technique. 

The  application  of  these  methods  to  real  as  opposed  to  simulated  data  introduced  elements 
not  considered  in  References  2  to  4.  One  of  these  elements  is  the  location  of  the  accelerometers 
relative  to  the  aircraft  centre  of  gravity.  The  location  is  rarely  coincident  with  the  centre  of 
gravity  and  this  introduces  certain  contributions  to  the  accelerometer  readings  due  to  aircraft 
rotational  rates  and  accelerations.  These  can  be  allowed  for  provided  the  position  of  the 
accelerometer  relative  to  the  centre  of  gravity  is  known.  Alternatively,  the  kinematic  equations 
can  be  extended  to  include  these  offsets  as  unknowns  to  be  estimated.  The  latter  approach  is 
developed  further  in  this  note.  The  present  study  also  extends  the  earlier  results  by  inclusion 
of  calibration  slope  errors  in  addition  to  the  bias  errors  treated  in  References  2  to  4.  This  exten¬ 
sion  as  well  as  the  effects  of  centre  of  gravity  offset  errors  arc  initially  studied  using  simulated 
data. 

Another  element  to  be  considered  with  real  data  is  the  presence  of  dynamic  lags  in  the 
measurement  channels.  Lags  may.  for  example,  be  introduced  by  the  signal  conditioning  elec¬ 
tronics.  Unless  this  is  identical  for  each  channel  it  cannot  be  assumed  that  the  measurements 
will  be  precisely  in  phase.  Phase  lags  of  up  to  around  0- 1  second  may  need  to  be  accounted  for. 
Another  source  of  lag  is  frequency  response  limitations  in  the  transducer  itself  or  in  the  measure¬ 
ment  system.  This  is  particularly  so  with  the  air  data  system  (e.g.  airspeed,  altitude)  which 
requires  movement  of  air  in  and  out  of  reservoirs  through  lubes  of  different  lengths.  Lags  intro¬ 
duced  in  this  way  may  amount  to  several  tenths  of  a  second  or  more  and  their  influence  on  the 
match  between  calculated  and  measured  time  histories  and.  consequently,  on  estimated  bia.s 
and  scale  factors  may  be  considerable.  Provided  measurement  system  frequency  response  band¬ 
width  is  sufficiently  wide  (typically  natural  frequency  S  to  10  times  the  frequency  being  measured), 
lags  can  be  approximated  by  simple  time  shifts.  L^  (or  time  shift)  effects  are  studied  and  quanti¬ 
fied  in  this  note  and  a  strategy  developed,  using  simulated  data,  to  determine  lag  values.  Some 
discussion  of  the  effects  of  lags  in  airspeed  measurement,  and  also  of  accelerometer  offsets, 
on  estimated  parameters  can  also  be  found  in  Reference  6. 

In  Section  2  of  this  report  the  compatibility  checking  method  is  reviewed  briefly  and  exten¬ 
sions  and  modifications  relevant  to  the  current  investigation  arc  described.  Results  with  these 
modifications,  including  centre  of  gravity  location  and  measurement  lag  determination  using 


simulated  data,  are  presented  in  the  next  Section.  Finally,  the  techniques  developed  are  applied 
to  flight  data. 


2.  THEORETICAL  BACKGROUND 

The  aircran  kinematic  equations  of  motion,  modified  to  account  for  centre  of  gravity 
offsets,  are  summarized  in  this  Section.  This  is  followed  by  a  brief  theoretical  discussion  of  the 
Maximum  Likelihood  (ML)  method  for  solution  of  the  compatibility  checking  problem  and  a 
more  detailed  treatment  of  the  sensitivity  matrix. 


2.1  Aircraft  Kiaematk  Equations 

The  full  set  of  non-linear  equations  relating  the  position,  velocity  and  acceleration  of  an 
aircraft  with  reference  to  a  set  of  flat  earth  axes  is  summarized  in  Reference  3.  In  the  present 
study  only  the  three  degrees  of  freedom  set,  representing  longitudinal  motions,  is  of  interest. 
The  system  under  consideration  may  be  summarized  as  follows: 


(i)  State  vector,  x  =  {u,  w,  A]''’ 

The  state  equations  written  in  body  axes  (x,  z)  fixed  in  the  aircraft,  with  origin  at  the  centre 
of  gravity  can  be  considered  to  be  exact  (Ref  7). 

(}  =  ~qw  +  axcii—g  sin  0  ' 

M'  =■  qu  +  azrr.  +g  cos  t> 

—  u  sin  0  —  w  cos  ti 

where 

u.  H'  are  linear  velocities  in  the  x,  z  directions, 

A,  q  are  pitch  attitude  and  angular  rate. 

azcii  are  accelerations  at  the  centre  of  gravity  in  the  x,  : 
directions,  and 

g  is  the  gravitational  constant. 


Accelerometers  are  seldom  located  exactly  at  the  aircraft  centre  of  gravity.  If  the  accelero¬ 
meter  pack  is  assumed  to  be  positioned  in  the  plane  of  symmetry  but  offMt  a  distance  xm  in  the 
-Y-direction  and  z>  in  the  z-diiection,  then  the  relationship  between  the  accelerations  at  the  r.g.  and 
those  indicated  at  the  accelerometer  location,  ax,  and  az,.  can  be  written  as  follows: 


axi  i;  ^  (ixn+q^x,  qz» 
azri!  —  azi,+qxa+q*ZM 


(2) 


These  non-linear  relations  require  a  knowledge  of  x*  and  z«  for  evaluation.  Alternatively. 
.Y(  and  Zt  may  be  treated  as  unknowns  and  added  to  the  list  of  parameters  to  be  identified.  The 
latter  approach  is  developed  further  in  this  report.  In  addition,  evaluation  of  Equation  (2)  requires 
knowledge  of  the  pitch  angular  acceleration,  4.  This  is  not  normally  one  of  the  quantities  measured 
and  the  present  approach  is  to  differentiate  records  of  pitch  rate,  q,  using  moving  least  squares 
smoothing  (Ref  8)  in  order  to  avoid  the  introduction  of  phase  shifts.  This  approach  gives  good 
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results  provided  care  is  taken  to  optimize  the  order  of  fit  and  the  number  of  points  used  in  the 
fit,  and  provided  also  that  noise  levels  on  q  are  not  too  high,  a  condition  which  is  fortunately 
satisfied  due  to  the  quality  of  pitch  rate  gyro  instrumentation  used. 


(li)  Input  vector,  u  =  [oxco,  ozco, 

In  the  compatibility  checking  approach,  values  of  axca,  ozco  t^nd  q  in  Equation  (t)  are 
treated  as  inputs.  While  Equation  (2)  enables  corrections  to  be  made  for  accelerometer  offset 
relative  to  the  centre  of  gravity,  measured  quantities  are  also  subject  to  errors  due  to  calibration 
errors,  instrument  bias  errors  and  random  noise.  These  are  modelled  as  follows: 


axt  =  (1 +  ^x)oxai+b»x+n»x 
=  (I  +A(x)aZiii  +  Iixx'('nxx 
^  =  (I  +  ^q)^in  + Aq  +  riq 


(3) 


where 


Aqx,  Au  and  Aq  represent  scale  factor  errors, 
btx,  and  Aq  represent  instrument  bias  errors, 

"ax,  ntz  and  rtq  represent  random  noise  and 
subscript  m  refers  to  measured  values. 

When  Equations  (2)  and  (3)  are  substituted  into  Equation  (1)  the  scale  factor  and  bias  errors 
as  well  as  the  accelerometer  positions  appear  explicitly  as  unknown  parameters  while  the  random 
noise  components  give  rise  to  process  noise.  Since  accelerometers  and  gyros  are  usually  the  most 
accurate  instruments  used  in  flight  testing,  the  noise  levels  are  in  general  relatively  small  and 
consequently  are  neglected  in  the  present  treatment.  This,  together  with  the  fact  that  the  state 
equations  can  be  considered  to  be  exact,  ensures  that  the  state  equations  are  free  of  process  noise. 


(Hi)  Output  vector,  z  =  [V,  «,  0,  AF 

Measurements  of  the  output  quantities,  viz.  airspeed,  K,  angle  of  attack,  a,  pitch  attitude,  0 
and  altitude,  h,  are  corrupted  by  scale  factor  errors,  biases  and  random  noise.  The  equations 
used  for  calculating  the  outputs  are  as  follows: 

►'out  =  (l+AvKM*  +  W*)>'*  +  Av  V 

«oa<  =  ( I  +  A,)tan--»l(»’-^x.)/«l + A,  1  (4) 

®oot  =  (I  +  A#)tf+A#  j 

Aou«  =  (I  +  Ah)A+Aji  / 


The  term  qxju  in  the  a  equation  is  a  correction  due  to  known  sensor  position,  x,,  ahead 
of  the  centre  of  gravity. 

The  biases.  A,  and  scale  factors.  A,  in  equation  (4)  are  added  to  the  list  of  unknown  para¬ 
meters  introduced  previously.  Further,  the  initial  conditions  of  the  state  equations  (I)  are  also 
unknown  and  complete  the  list  of  parameters  to  be  determined.  The  list  can  be  reduced  slightly 
by  noting  that,  since  no  absolute  reference  for  height.  A,  is  implied  in  equation  (I),  the  initial 
height  A(o),  and  bias  As  can  be  assumed  to  be  zero.  The  final  unknown  parameter  vector  is: 

(  *=  IA»x,  Asx*  A«f,  Am>  Aq,  Aq,  A»,  by,  A„  A„  Aq,  Aq,  Aq,  m(o),  w(o),  A(e),  x»,  2«F  (5) 
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2.2  Ilnory 

The  system  described  in  the  previous  section  form  a  set  of  non-linear  dynamic  equations 
of  the  form: 


x(r)  =  f(x(r),  o(/),  0 

(6) 

*(0  =  g(x(<).  B(0.  f) +■(»■) 

(7) 

The  measurements,  z(j),  are  made  of  a  finite  number  of  time  points,  u,  and  are  subject  to 
random  measurement  noise,  n,  in  addition  to  the  systematic  errors  represented  through  the 
parameter  vector,  f 

Since  there  is  no  process  noise,  an  estimate  of  the  state,  x,  is  obtained  by  integration  of 
equation  (6)  given  the  estimated  value  of  the  parameter  vector.  Parameter  estimation  is  accom¬ 
plished  using  a  Maximum  Likelihood  (ML)  method.  Assuming  that  the  measurement  noise  is 
zero  mean  uncorrelated  and  Gaussian,  this  is  equivalent  to  the  minimization  of  the  cost 
functional  (Ref.  9). 


N 

J(i,  R)  =  t  A  log  1  R  1  +  i  y  vT(i)  R-i  v(i) 


(8) 


where  v(i)  =  z(/)— g(*(i),  u(/),  =  z(i)— i(f,  i)  is  the  vector  of  residuals,  and  R  is  a  weighting 

covariance  matrix. 

An  estimate  for  R  is  first  obtained  by  minimizing  J  (equation(8))  with  respect  to  R,  giving 

N 

*  =  (9) 

/•=  1 

An  updated  estimate  for  f  then  follows  by  minimizing  /  with  respect  to  (.  This  revised 
estimate  is  used  for  a  new  estimate  of  R  and  the  iteration  procedure  continued  until  convergence. 
The  final  estimates  are  asymptotically  consistent  and  efiSdent  with  covariance  of  the  estimates 
approaching  the  Cramer-Rao  bound  given  by: 


(10) 


where  A  is  the  sensitivity  matrix  whose  elements  are  the  partial  derivatives  of  the  elements  of 
the  estimated  output  vector,  2(f,  i)  with  respect  to  elements  of  the  parameter  vector,  (. 

A(/)  =  Vj(l(f,/))  (II) 

Further  details  are  provided  in  Reference  3. 


covariance  (f)  = 


N 

1 

U=  1 


AT(i)R-iA(f) 


2.3  Scildvity  Matrix 

A  critical  part  of  the  ML  procedure,  both  for  the  parameter  estimation  and  for  the 
calculation  of  the  covariaricc  of  the  estimates,  is  the  evaluation  of  the  sensitivity  matrix.  For 
the  current  problem  the  output  vector  (equation  (4))  is  of  dimension  4  and  the  parameter  vector 
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(equation  (5))  is  of  dimension  18  resulting  in  a  sensitivity  matrix  of  dimension  4  by  18.  The  (/,  j) 
element  is  the  partial  derivative  of  the  tth  component  of  the  output  vector  with  respect  to  the 
/th  component  of  the  parameter  vector.  The  complete  sensitivity  matrix.  A,  for  the  present 
system  is  summarized  in  Appendix  I  and  formulae  for  the  evaluation  of  the  various  elements 
given. 

in  order  to  evaluate  the  elements  of  A  as  functions  of  time,  it  can  be  seen  that  not  only 
the  state  u.  w,  0,  h  but  also  partial  derivatives  of  the  state  are  required  a.s  functions 

of  time,  where  represents  a  particular  element  of  the  parameter  vector.  Equations  for  the 
required  partial  derivatives  follow  from  the  state  equations  by  first  differentiating  with  respect 
to  and  then  reversing  the  order  of  the  partial  derivative  and  the  time  derivative.  The  total 
number  of  equations,  including  the  state  equations,  which  need  to  be  integrated  simultaneously 
amounts  to  39  and  these  are  listed  in  Appendix  2  together  with  the  relevant  initial  conditions. 

For  computing  purposes  derivative  values  for  all  variables  to  be  integrated  are  evaluated  in 
subroutine  DERIVS  while  manipulations  required  for  calculation  of  A  are  performed  in  sub¬ 
routine  SENS.  Further  details  are  given  in  Reference  3.  The  main  program  COMPAT  is  written 
so  that  the  user  may  choose  any  desired  subset  of  the  unknown  parameters  for  extraction.  The 
program  recognizes  the  required  subset  from  the  input  file  and  makes  the  necessary  adjustments 
automatically. 


3.  RESULTS  WITH  SIMULATED  DATA 

In  this  Section  the  method  of  the  previous  Section  is  used  to  study  those  practical  aspects 
described  in  the  introduction.  These  include  in  turn,  the  addition  of  scale  factor  errors,  accelero¬ 
meter  offset  relative  to  the  centre  of  gravity  and  the  effects  of  lags  in  the  measurements.  Finally, 
a  strategy  is  developed  for  the  determination  of  lags  and  tested  using  simulated  time  histories. 


3.1  Manoeuvre  Description 

Two  manoeuvre  shapes  were  used  in  the  course  of  this  study,  both  being  representative  of 
rollercoaster  type  of  motion  in  the  longitudinal  plane,  i.e.  sequences  of  pull-ups  and  push-overs: 


Manoeuvre  /  (Fig.  I). — This  is  similar  to  the  manoeuvre  used  in  Reference  4.  Analytical 
expressions  for  the  inputs  ax,  az  and  q  are  chosen  to  simulate  a  real  manoeuvre.  The  higher 
frequency  oscillations,  particularly  in  the  pitch  rate.  q.  time  history,  represents  the  short  period 
oscillation.  Maximum  deviations  are  seen  to  be  about  :l  0-5  m/s  (0-05  g)  for  ax. 
+  I0m/s  (Ig)  for  az  and  ±0-2rad/s  (lldeg/$)  for  q.  The  corresponding  outputs  (noise  free) 
used  in  the  current  study,  namely  F.  «  and  0  arc  also  shown  in  Figure  1 . 


Manoeuvre  2  (Fig.  2). — The  second  manoeuvre  used  is  a  close  representation  of  a  real 
manoeuvre.  Actual  measured  elevator  time  histories  were  used  as  input  to  a  mathematical 
model  of  the  aircraft  in  question  and  the  dynamic  response  calculations  of  ax,  az  and  q  produced 
by  the  model  then  used  as  inputs  for  the  current  exercise.  This  manoeuvre,  although  more 
realistic  than  manoeuvre  I,  produces  somewhat  smaller  deviations  in  accelerations  and  pitch 
rate. 
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An  auxiliary  program  was  used  to  add  the  effects  of  bias  and  scale  factor  errors  ot'  measure¬ 
ments  of  inputs  and  outputs  and  of  random  noise  on  output  measurements.  Three  different 
levels  of  noise  have  been  used  at  one  time  or  another  and  these  are  presented  in  Table  I  below. 


TABLE  I 

Noise  Standard  Errors  for  Simulated  Data 


Standard  Error 

Level  1 

Level  2 

Level  3 

<7(  F),  m/s 

1  0 

0  1 

0  5 

aix),  rad 

0  002 

0  001 

0  002 

aib),  rad 

001 

0  001 

0  002 

Level  I  represents  a  high  noise  case,  level  2  low  noise  and  level  3  an  intermediate  case. 
Measured  outputs  shown  in  Figure  2  represent  the  level  I  noise  case.  Also  shown  in  the  figure 
is  a  time  history  of  pitch  acceleration,  q.  This  has  been  derived  from  the  pitch  rate,  q,  history 
by  numerical  differentiation  as  explained  in  Section  2.1.  The  maximum  pitch  rate  apparent  in 
Figure  2  is  approximately  0-08  rad/s  while  maximum  pitch  acceleration  is  about  0  ISrad/s‘^. 
For  centre  of  gravity  offsets,  .v,  and  2a,  of  say  I  metre  these  correspond  to  corrections  to 
accelerometer  readings  (equation  (2))  of  0  006  m/s*  and  0- 15  m/s*  respectively.  Clearly  the  q 
contributions  are  dominant  and  can  make  significant  contributions  to  the  u  and  w  state  equations, 
comparable  at  times  to  the  contributions  of  the  bias  terms,  bax  and  bai.  which  are  set  at  0  - 1  m/s* 
in  the  simulated  runs  of  the  following  Section. 


3.2  Results  with  Scale  Factor  Errors  and  Accelerometer  Offsets 

All  simulation  results  reported  in  this  and  the  next  Section  were  obtained  using  40  seconds 
of  record  at  40  samples  per  second,  shown  to  give  satisfactory  results  in  References  3  and  4. 
CPU  time  on  the  DEC-System  10  computer  for  10  iterations  of  the  minimizing  procedure  was 
typically  20  minutes,  varying  only  slightly  with  the  number  of  parameters  extracted.  Only  the 
three  outputs  F,  <x  and  f>  were  matched  since  Reference  3  suggests  that  relatively  minor  improve¬ 
ments  are  obtained  by  the  addition  of  height  h,  as  an  output. 

The  results  including  scale  factor  and  accelerometer  offset  extraction  are  summarized  in 
Tables  2  and  3  .  The  full  parameter  vector  considered  here  is: 

i  -  (b«x.  bait,  A,|,  b(|.  Av,  by,  Aj,,  b_,,  b/f,  ufo),  m'(u),  .v».  2x] 

and  includes  the  scale  factors  A,,,  Ay,  A^  and  accelerometer  offsets  ATa.  z%.  Normally  scale  factor 
errors  for  gyros  and  accelerometers  can  be  expected  to  be  small,  although  in  the  current  exercise 
A,,  has  been  retained  due  to  inconsistencies  noted  in  the  pitch  rate  gyro  calibrations  for  the  flight- 
data  to  be  considered  in  Section  4.  A  priori  estimates  for  all  parameters  except  the  initial  con¬ 
ditions  were  set  to  zero.  Initial  condition  estimates  were  derived  from  the  simulated  outputs. 

The  results  in  Table  2  were  extracted  using  manoeuvre  2  and  low  noise  levels.  In  the  first 
column  no  accelerometer  offsets  are  present  but  all  other  parameters  are  included.  Most  para¬ 
meters  are  well  identified  and  the  Cramer-Rao  bounds  (in  brackets)  are  good  indicators  of  their 
accuracy  as  noted  in  References  3  and  4.  It  is  particularly  pleasing  that  the  scale  factors  Aq,  Ay 
and  A^  are  very  precisely  identified  and  present  no  special  problems.  The  second  column  includes 
accelerometer  offsets  of  ±  I  m  in  both  the  x  and  z  directions.  Their  inclusion  leads  to  no  notice¬ 
able  deterioration  in  the  accuracy  of  the  other  parameters.  However,  while  and  2.  arc 
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reasonably  well  identified  the  accuracy,  particularly  of  as  indicated  by  the  Cramer-Rao 
bound,  is  relatively  poor.  This  is  most  probably  associated  with  the  small  overall  contributions 
the  q  terms  make  in  the  equations  of  motion  for  this  manoeuvre,  compared  to  the  other  bias  and 
scale  factor  errors.  The  absence  of  other  errors  can  lead  to  an  improvement  in  accuracy  for 
JTa  and  Za  as  shown  in  column  3  of  Table  2. 


TABLE  2 

Results  with  Scale  E^rrors  and  C.G.  Offsets 


Parameter 

True  value 

1 

2 

3 

fiax,  m/s® 

01 

0-114  (0-003) 

0-095  (0-003) 

0-098  (0-003) 

baz,  m/s® 

01 

0-102  (0-0007) 

0-100  (0-0007) 

— 

Aq 

001 

0-0098  (0-0002) 

0-0101  (0-0002) 

— 

bq,  rad/s 

0  002 

0-0020(0-000002) 

0-0020(0-000002) 

— 

Ay 

01 

0-100  (0-0003) 

0-099  (0-0004) 

— 

m/s 

10 

0-82  (0-07) 

1-16  (0-08) 

1-02  (0-03) 

-'a 

01 

0-099  (0-0005) 

0-100  (0-0007) 

— 

b^,  rad 

0  002 

0-0022(0-0002) 

0-0020(0-0002) 

— 

be,  rad 

001 

0-0086(0-0003) 

0-0105  (0-0003) 

0-0102(0-0003) 

u{o),  m/s 

193-76 

193-89  (0-04) 

193-74  (0-05) 

193-74  (0-03) 

tv(o),  m/s 

28-85 

28-88  (0-02) 

28-85  (0-02) 

28-85  (0-01) 

6{o),  rad 

0-148 

0-148  (0-003) 

0-147  (0-0003) 

0-148  (0-0003) 

x»,  m 

-1-0 

— 

-0-90  (0-20) 

-1-04  (0-13) 

z»,  m 

l-O 

0  97  (0-09) 

1-05  (0  07) 

Manoeuvre  2;  Noise  Level  2. 


Alternatively,  the  accuracy  of  extraction  of  Xa  and  Za  can  be  improved  by  increasing  the 
contributions  of  the  q  terms.  There  is  considerable  scope  for  g  this  since  the  manoeuvres 
studied  here  are  relatively  mild.  Thus  a  manoeuvre  with  large  rate  accelerations  could  provide  a 
method  for  obtaining  aircraft  centre  of  gravity  position  directly  from  flight  tests.  This  requires 
further  investigation  but  is  not  pursued  here. 

Table  3  presents  some  results  showing  the  effects  of  bias  and  scale  factor  magnitudes  and 
random  noise  levels  on  the  accuracy  of  extraction  of  the  various  parameters.  The  first  two  columns 
both  have  level  2  (low)  noise  but  bias  and  scale  factor  errors  are  twice  as  large  in  the  first  column 
as  in  the  second.  The  Cramer-Rao  bounds,  determined  principally  by  the  noise  level,  arc  iden¬ 
tical  in  both  cases  and  most  parameters  are  identified  equally  well.  A  possible  exception  is  6v 
and  more  particularly  Xa,  accelerometer  offset  in  the  x-direction.  Both  of  these  appear  to  be 
biased  in  the  first  column  in  the  sense  that  their  deviations  from  the  true  values  are  well  outside 
the  predicted  Cramer-Rao  error  bounds. 

Columns  2,  3  and  4  of  Table  3  compare  results  with  low,  high  and  intermediate  noise  levels 
respectively,  but  with  the  lower  bias  and  scale  factor  errors.  The  deterioration  in  accuracy, 
in  particular  the  Cramer-Rao  bound,  as  noise  levels  increase  is  obvious.  In  column  4  the  results 
arc  the  mean  of  eight  separate  runs  and  the  standard  deviations  are  closely  approximated  by  the 
Cramer-Rao  bounds. 

In  summary,  this  Section  has  shown  that  the  presence  of  scale  factors  does  not  present  any 
difficulties  to  the  compatibility  checking  algorithm,  which  is  capable  of  calculating  their  values 
with  very  good  accuracy.  For  the  relatively  mild  manoeuvre  considered,  the  accelerometer 
offsets  were  only  identified  with  moderate  accuracy,  however  the  presence  of  offsets  did  not 
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TABLE  3 


have  any  adverse  effects  on  the  other  parameters.  Finally,  increases  in  noise  levels  lead  to 
significant  degradation  in  accuracy  but  increased  bias  and  scale  factor  magnitudes  do  not  in 
general  affect  the  relative  accuracy  of  the  identified  values. 


33  Results  with  Lags 

The  questions  addressed  in  this  Section  include  what  the  effects  of  lags  arc  on  the  identified 
parameter  values  and  how  measurement  lags  affect  the  fit  errors  between  calculated  and  measured 
time  histories.  In  order  to  answer  these  questions  a  systematic  study  was  carried  out  whereby 
a  lag  of  ±  0  05  s  was  imposed  on  each  of  the  inputs  (a.T.  a:,  q)  and  outputs  (F,  *.  <l)  in  turn. 
For  this  study  level  2  (low)  noise  level  was  used  and  the  parameter  vector  was; 

f  =  [btx,  bxz,  Ag.  bg.  6v.  b„  b/»,  mo),  nio). 

Both  manoeuvres  I  and  2  were  used  and,  as  may  be  expected,  the  main  parameters  to  be 
affected  were  those  that  were  found  to  be  least  accurately  determined  in  the  previous  Section 
(also  see  Ref.  3),  viz.  b^x,  h*.  b^  and  The  results  arc  presented  qualitatively  in  Table  4  below, 
where  a  major  change  (100%  or  more)  in  a  parameter  is  denoted  by  a  large  x  while  a  relatively 


TABLE  4 

Effects  of  Lags  on  Identified  Parametets 


N^arameter 

Lag  inS. 

b%x 

by 

bx 

b» 

ax 

az 

X 

X 

9 

X 

X 

X 

X 

V 

+ 

X 

S 

+ 

X 

4^ 

9 

+ 

+ 

Manoeuvre  1,2;  Noise  Level  2. 


minor  change  is  indicated  by  a  small  + .  Thus  a  lag  of  0  05  second  in  az  leads  to  major  changes 
in  bxx  and  bf  and  a  small  but  significant  change  in  bx.  Note  that  there  are  only  two  major  changes 
in  each  column  and  they  are  in  general  of  opposite  sign.  For  example  is  strongly  affected 
by  a  lag  in  az  and  similarly  affected,  in  the  opposite  direction,  by  a  lag  in  q.  This  suggests  that 
the  main  influence  of  b^x  is  the  relative  phase  between  az  and  q.  Similarly  for  bf  while  for  bv 
the  phase  between  q  and  a  is  important  and  for  b„  that  between  q  and  K  TTius  a  lag  in  q  affects 
all  four  parameters  in  a  major  way,  a  lag  in  az  affects  b^x  and  b,.  while  lags  in  V  and  *  influence 
only  b,  and  by  respectively.  A  lag  in  a.x  has  no  noticeable  effect  on  any  parameter  and  a  lag  in 
n  affects  b»x,  by  and  b,  in  a  minor  way. 
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Results  showing  the  influence  of  lags  on  the  V,  a  and  6  fit  errors  are  summarized  in  Table  5. 
The  numbers  relate  to  manoeuvre  2  and  are  weighted  RMS  fit  errors  obtained  after  5  iterations 
of  the  minimization  algorithm.  For  no  lags,  values  dose  to  1  would  be  obtained  in  all  cases. 
Small  variations  may  be  expected  due  to  differences  in  the  random  noise  sequences  used  in  each 
case.  The  results  show  clearly  that  the  fit  error  for  V  is  affected  only  by  the  relative  phase  of 
q  and  V  while  the  P-fit  depends  on  the  relative  phase  of  q  and  6.  The  «-fit  on  the  other  hand  is 
influenced  by  the  relative  phases  of  q,  a:  and  o.  Once  again  the  q  lag  appears  to  be  dominant 
while  ax  lag  has  no  discernible  effect  on  the  fit  errors.  On  the  other  hand  8  lag  with  respect  to  q, 


TABLE  5 

Fit  Errors  Obtained  After  5  Iterations 


O  OS  s  Lag  in 

F-fit 

a-fit 

9.fit 

ax 

103 

1  01 

I  01 

az 

1-04 

1-47 

1  04 

4 

1-33 

2-41 

4-30 

F 

1  25 

0  98 

101 

a 

0  96 

1  94 

1  03 

9 

I -01 

0  93 

4  48 

Manoeuvre  2;  Noise  Level  2. 


while  influencing  some  parameters  in  a  minor  way  (Table  4),  has  a  large  effect  on  the  9-fit.  At 
the  same  time  V  lag  w.r.t.  q  has  a  significant  influence  on  the  F-fit  corresponding  with  its  effect 
on  in  Table  4.  Finally,  Table  4  suggests  that  the  relative  phase  of  az  and  «  is  relatively  un¬ 
important  so  that  the  a-fit  results  in  Table  5  can  be  assumed  to  reflect  q—az  and  q—a  phase 
differences  independently,  the  a—az  phase  relationship  being  of  secondary  importance.  These 
observations  lead  to  the  following  strategy  for  determining  relative  phases: 

1 .  Take  q  as  reference  signal. 

2.  Align  6  with  q  by  minimizing  9-fit  w.r.t.  9  lag. 

3.  Align  a  with  q  by  minimizing  a-fit  w.r.t.  a  lag. 

4.  Align  az  with  q  by  minimizing  a-fit  w.r.t.  az  lag. 

5.  Align  y  with  q  by  minimizing  F-fil  w.r.t.  F  lag. 

The  above  steps  2  to  S  can  be  iterated  to  check  possible  coupling  effects.  In  particular 
repetition  of  steps  3  and  4  can  be  repeated  to  confirm  that  az  and  a  can  be  independently  aligned 
with  q  and  that  the  relative  alignment  of  az—»  is  unimportant. 

In  order  to  check  the  above  strategy  a  set  of  manoeuvre  1  time  histories  with  level  1  (high) 
noise  was  generated  and  then  shifted  relative  to  one  another  to  produce  the  relative  lags  shown 
in  Table  6.  The  parameter  vector  used  in  this  simulation  included  the  scale  factors  Av  and  A,, 


10 


TABLE  t 


RdalhcLaci 


Lag  in 

Seconds 

Sample  intervals 

ax 

005 

2 

02 

0  05 

2 

8 

0 

0 

F 

0-275 

II 

at 

-0  1 

-4 

8 

-0-05 

-2 

but  not  Ao,  together  with  biaiet  and  initial  conditions  as  before.  Starting  with  step  2  of  the 
strategy,  the  8  record  was  qrstematkally  shifted  and  each  time  the  compatibility  checking  pro¬ 
gram  used  to  match  records.  The  shift  which  produced  the  least  0-fit  error  after  five  iterations 
of  the  algorithm  established  the  relative  lag  of  the  9  record  with  respect  to  the  ^  reference.  This 
same  general  procedure  was  repeated  for  the  s-record  (step  3)  and  so  on.  At  each  stage  the 
appropriate  time  history  was  shifted  by  the  identified  number  of  lag  time  intervals  before  pro¬ 
ceeding  to  the  next  step. 

The  results  of  the  above  procedure  are  summarized  in  Figures  3  and  4.  In  Figure  3  the 
minimum  point  of  the  9-fit  error  identifies  the  correct  8  lag  value  (-2)  and  similarly  the  s-fit 
error  identifies  the  a  lag  value  (—4).  Also  shown  are  the  variations  in  the  other  fit  errors  as 
9  lag  or  a  lag  are  varied.  As  expected,  these  variations  are  small  by  comparison.  In  Figure  4 
(upper)  the  correct  az  lag  value  (2)  is  obtained  by  minimizing  the  «-fit  error  with  respect  to 
or  lag.  It  can  be  assumed  that  <ur  lag  is  the  same  as  az  lag  although,  as  noted  earlier,  the  value 
of  ax  lag  does  not  significantly  affect  the  results.  Finally,  Figure  4  (lower)  determines  F  lag  ( 1 1 ) 
by  minimizing  the  F-fit  error  with  respect  to  F  lag.  Both  minima  in  Figure  4  are  weaker  than 
those  shown  in  Figure  3.  but  the  lag  values  are,  nevertheless,  unambiguously  identified.  With 
the  determination  of  F  lag  the  process  is  complete  and  examination  of  Figures  3  and  4  shows 
that  the  final  correct  set  of  lags  corresponds  to  the  absolute  minimum  lit  errors  for  each  of  the 
9,  a  and  F  fits.  Further,  it  has  not  been  found  necessary  to  iterate  the  procedure. 

These  results  with  simulated  data  show  that  it  is  possible  to  determine  lag  values  accurately 
by  systematic  shifting  of  input  and  output  time  histories  using  the  q  history  as  reference. 
Quantitative  differences  may  occur  with  different  manoeuvre  shapes  but  the  present  results 
show  good  qualitative  agreement  in  results  from  manoeuvre  I  (Figures  3  and  4)  and 
manoeuvre  2  (Table  S). 


4.  LAG  RESULTS  WITH  FLIGHT  DATA 

The  practical  considerations  studied  in  the  previous  Section  are  applied  here  to  data  obtained, 
on  an  opportunity  basis,  during  performance  flight  triah  of  a  delta  wing  fighter  aircraft  by  the 
RAAF  Aircraft  Research  and  Devefopment  Unit  (ARDU).  The  manoeuvre  and  data  acquired 
will  first  be  described  followed  by  diacusHon  of  some  of  the  steps  required  in  preparation  for 
application  of  the  compatibility  checking  program.  The  results  for  instrumentation  lag  will 
then  be  presented. 


II 


4.1  M—ofyrr  mi  P«U  Ae^iiiiriM 

The  manoeuvre  under  consideration  is  a  roller  coaster  manoeuvre  commencing  at  an  altitude 
or  approximately  33QOOft  and  a  Mach  number  of  0*6S.  Lateral  motions  were  minimal  so 
that  the  manoeuvre  could  be  considered  to  be  purely  longitudinal.  Hight  test  instrumentation 
included  a  pitch  rate  gyro,  pitch  attitude  gyro  and  accelerometers  to  measure  nomtal  and  axial 
accelerations.  A  nose-probe  fitted  with  a  hemisphericai,  S-port  flow  angle  sensor  was  used  for 
angle  of  attack  measurements.  This  design  placed  pressure  transducers  within  the  probe  and  as 
close  to  the  probe  head  as  possible,  thereby  reducing  time  lags  due  to  the  tubing  to  negligible 
amounts.  The  probe  was  calibrated  in  the  ARL  transonic  wind  tunnel.  Airspeed  measurements 
were  only  available  through  the  aircraft  air  data  system  and  presumably  subject  to  significant 
time  lags.  Finally,  pressure  altitude,  although  not  used  in  determining  lags,  was  also  available. 

Data  were  digitally  recorded  (12  bit  accuracy)  on  magnetic  tape  via  the  Aircraft  Flight  Test 
Recording  and  Analysis  System  (AFTRAS — R^.  10)  used  by  ARDU,  at  a  sampling  rate  of 
60  per  second.  Instrument  lags  and  lags  introduced  by  the  data  conditioning  prior  to  recording 
were  unknown.  Copies,  on  seven  track  tape,  of  the  ARDU  flight  tapes  were  provided  to  ARL 
and  processed  as  described  in  Reference  1 1  to  obtain  the  required  time  histories  in  Engineering 
units.  Time  histories  for  the  manoeuvre  under  consideration  are  shown  in  Figure  S.  The  roller 
coaster  motion  during  the  40  seconds  of  flight  shown  is  apparent,  with  an  initial  decrease  in  air¬ 
speed  as  the  aircraft  climbs  followed  by  an  increase  as  the  aircraft  descends  and  so  on.  Pitch 
dampers  were  not  operating  during  this  flight  as  evidenced  by  the  short  period  oscillations 
particularly  noticeable  in  the  az,  q  and  a  records. 

From  Figure  5  it  can  be  seen  that  the  range  of  values  of  ax  is  relatively  small  spanning  about 
±  0-3  m/s*  (0  03  g).  Meanwhile  the  range  of  the  az  measurements  is  around  ±  10  m/s*  (I  g) 
and  q  spans  approximately  ±0-1 1  rad/s  (6 deg/s).  These  values  represent  a  relatively  gentle 
manoeuvre.  Noise  of  the  q  and  az  measurements  is  not  discernible  on  the  traces  while  that  on 
ax  appears  to  be  confined  to  a  range  of  ±  O-OS  m/s*  (O-OOSg).  Larger  scale  plots  suggest  a 
similar  noise  level  (O-OSm/s*)  on  az  and  ±0-002rad/s (0-1  deg/s)  on  q.  Noise  is  much  more 
obvious  on  the  y  and  fi  outputs  being  about  ±  I  m/s  (2  knots)  and  ±  25  f)  (8  m)  respectively. 
However,  for  a  and  B,  noise  is  difficult  to  estimate  on  the  scales  given,  but  from  larger  scale  plots 
random  noise  was  estimated  to  be  at  most  ±  0-002  rad  (0  1°)  for  the  a  and  B  records.  Noise 
bandwidth  was  unknown  in  all  cases. 


4.2  Preliaiinary  CaIcnfartioM 

The  pitch  rate  acceleration,  4.  required  in  order  to  take  into  account  accelerometer  offset 
from  centre  of  gravity  was  calculated  by  numerical  differentiation  using  a  moving  least  squares 
fit  as  explained  in  Section  2. 1 .  Some  experimentation  was  necessary  with  the  number  of  points 
smoothed  and  the  order  of  the  polynomial  fit  requited  to  obtain  smoothly  varying  results. 
Figure  6  shows  the  original  pitch  rate  time  history,  q,  and  the  pitch  acceleration  record,  4,  as 
derived  using  a  fourth  order  polynomial  smoothing  40  points  at  a  time.  The  results  appear  to 
be  quite  satisfactory. 

A  priori  estimates  of  the  accelerometer  offsets  relative  to  the  c.g.  were  n»de  using  the  known 
position  of  the  accelerometers  and  calculated  values  of  the  c.g.  location.  The  latter  were  obtained 
using  a  mathematical  model  of  the  aircraf)  (Ref.  12),  given  the  aircraft  configuration,  initial 
weight  and  fuel  gone.  This  resulted  in  an  estimate  for  of  0  - 1 1  and  for  z*  of  0  -  S3  m. 

The  output  vector  in  the  compatibility  checking  program  itKiuded  F,  «  and  B  but  not  h. 
The  diagonal  weighting  nuitrix,  R,  in  equation  (8),  was  baaed  on  noise  standard  deviations  of 
I  m/s  for  y  and  0-001  rad  for  a  and  B.  ThcM  values  were  not  varied  during  the  S  iterations  of 
the  program  used  at  each  step  of  the  lag  dctermiruition  process.  The  results  have  been  found 
to  be  relatively  insensitive  to  the  value  of  the  R  matrix.  The  parameter  vector  used  reflects 
anticipated  systematic  errors  based  on  calibration  information  provided  by  ARDU.  In  particular 
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the  scale  factor  error  Aq  is  included  since  separate  calibrations  were  found  to  produce  significant 
differences  in  pitch  rate  gyro  calibration  slope.  The  complete  vector  was: 


f  —  lh»xt  fin.  Aq,  hq,  Av,  6v,  A„  h,,  Aq,  hq,  u{o),  wto),  6{o),  Zqp" 

It  was  noted  also  that  Xq  and  z*  were  not  important  in  the  current  manoeuvre  since  relatively 
small  pitch  rates  and  accelerations  were  present.  Consequently,  to  avoid  x.  and  Zq  taking  up 
unrealistic  values,  a  priori  weightings  were  specified  based  on  standard  deviations  of  0-S  and 
0*2  m  for  Xq  and  Zq  respectively.  All  other  parameters  were  allowed  to  vary  without  any  a  priori 
weighting.  In  some  cases  it  was  found  that  the  presence  of  a  priori  weights  on  parameters  could 
distort  the  result. 


4,3  Detenuinatiau  of  Lags 

The  results  using  the  procedure  outlined  in  Section  3.3  are  shown  in  Figures  7  and  8  and 
can  be  compared  with  the  corresponding  results  with  simulated  data  in  Figures  3  and  4.  Note 
that  the  lag  values  identified  in  Figures  7  and  8  must  be  compared  with  the  reference  q  lag  value, 
which,  for  practical  reasons,  is  6  rather  than  zero. 

As  with  the  simulated  data,  the  flight  data  (Fig.  7)  shows  the  8  and  a  fit  errors  to  be  strong 
functions  of  8  lag  and  a  lag  respectively.  The  8  lag  value  is  close  to  8  and  the  a  lag  value  is  3. 
Given  that  the  reference  q  lag  is  6,  this  means  that  the  «  signal  leads  the  q  signal  by  3  sample 
intervals  or  O  OS  seconds.  This  is  not  entirely  surprising  in  view  of  the  probe  design  as  discussed 
in  Section  4. 1 .  Oifferences  between  q  and  « lag  values  can  be  ascribed  to  transducer  characteristics 
and  signal  conditioning.  In  Figure  8  the  variations  of  a  and  F  fit  errors  with  az  lag  and  F  lag 
respectively  are  relatively  weak  but  nevertheless,  as  with  the  simulated  dau  of  Figure  4,  clear 
minima  can  be  discerned.  Thus  the  az  lag  value  is  about  12  while  the  F  lag  value  is  about  28. 
It  may  be  conjectured  that  a  more  violent  manoeuvre  leading  to  more  rapid  variations  in  aircraft 
state  would  result  in  a  more  pronounced  minimum  in  the  F  fit  curve.  Ihis  has  yet  to  be  tested. 
The  present  results  show  that  lags  can  be  identified  from  flight  data  in  a  straight  forward,  syste¬ 
matic  manner.  As  with  the  simulated  data,  iteration  of  the  procedure  was  not  needed. 

The  length  of  data  required  to  achieve  good  resulu  would  depend  on  the  information  present 
in  the  available  data,  i.e.  the  manoeuvre  shape.  In  the  present  instant,  increasing  the  data  length 
to  70  s  made  little  difference  to  the  identified  lags.  However  variations,  particularly  of  a  fit  error 
with  dz  lag,  were  somewhat  weaker.  This  suggests  that  shorter  data  lengths  would  be  satisfactory 
in  identifying  the  lags,  although  the  accuracy  of  the  parameters  extracted  would  deteriorate 
with  reduced  record  length  (Ref.  3).  The  accuracy  of  extraction  of  the  unknown  instrumentation 
parameters  and  their  dependence  on  measurement  lags  will  be  addressed  in  a  subsequent  note. 

Finally,  the  match  achieved  after  S  iterations  between  the  calculated  and  measured  histories 
of  F, «  and  8,  with  all  lags  adjusted  in  accordance  with  their  values  determined  as  above,  is  shown 
in  Figure  9.  For  a  and  8  the  calculated  and  measured  responses  are  almost  indistinguishable 
with  root  mean  square  fit  errors  of  0-0016  and  0-0022  radians  respectively.  The  error  in  the  F 
match  is  more  discernible  with  a  root  mean  square  value  of  0-61  m/s.  The  measured  F  record 
also  shows  hints  of  flatness  at  its  peak  and  valley,  indicating  possible  friction  or  deadband  effects 
in  the  air  data  computer.  Such  non-linear  effects  have  not  been  included  in  the  model  and  may 
possibly  influence  the  identified  lag  values.  Ideally,  the  ability  of  the  suggested  procedure  to 
identify  lag  values  correctly  should  be  verified  against  independent  experimentally  determined 
values. 


5.  CONCLUDING  REMARKS 

Several  matters  relating  to  the  application  of  instrument  compatibility  checking  techniques 
to  real  data  have  been  studied  in  this  report.  These  include  the  pretence  of  calibration  or  scale 
errors  in  addition  to  bias  errors,  accelerometer  offsets  relative  to  aircraft  centre  of  gravity,  and 


13 


relative  phare  or  time  ihifU  between  the  measuremente.  Only  symmetric  motions  of  the  aircraft 
were  considered. 

The  compatibility  checking  method  reformulates  the  aircraft  kinematic  equations  so  as  to 
make  use  of  the  redundant  information  available  in  the  measurements  to  extract  information 
about  the  instrument  systematic  errors  and  accelerometer  offsets.  A  Maximum  Likelihood 
program  is  used  to  solve  the  resulting  non-linear  estimation  problem.  Neglect  of  process  noise 
is  justified  because  of  the  exactness  of  the  kinematic  equations  and  the  low  noise  on  the  measured 
inputs. 

The  results  with  simulated  data  demonstrated  that  the  compatibility  checking  algorithm 
was  capable  of  calculating  scale  factor  errors  to  very  good  accuracy  even  at  high  noire  levels, 
although  some  general  degradation  of  accuracy  with  increased  measurement  noire  was  evident. 
On  the  other  hand,  for  the  relatively  mild  manoeuvres  considered,  accelerometer  offsets  were 
only  identified  to  moderate  accuracy.  However,  with  manoeuvres  designed  to  produce  large 
rate  accelerations  the  present  technique  offers  the  potential  of  in-flight  determination  of  centre 
of  gravity  axial  and  vertical  location  and  may  be  worth  investigating  further. 

Simulated  manoeuvres  also  demonstrated  that  relatively  small  phase  shifts  between  measured 
time  histories  can  have  major  effeett  on  the  values  of  extracted  parameters  and  also  on  the 
match  achieved  between  cahnilated  and  measured  dau.  The  latter  effect  was  used  to  devise  a 
systematic  procedure  to  identify  the  unknown  relative  phases.  The  procedure  was  applied 
successfully  to  simulated  records. 

The  lag  determination  procedure  was  also  applied  to  a  40-second  segment  of  flight  dau 
represenUtive  of  a  roller-coaster  manoeuvre.  The  compatibility  checking  program  in  this  case 
made  allowance  for  calibration  and  bias  errors  and  also  for  accelerometer  offset.  The  instru- 
menution  was  of  moderately  good  accuracy  and,  for  the  given  manoeuvre  shape  and  length, 
a  credible  set  of  lag  values  was  clearly  identified.  In  particular,  the  airspeed  measurement  lagged 
the  reference  pitch  rate  measurement  by  22  ti^  intervals  or  approximately  0-37  seconds.  A 
check  against  an  independent,  experimenully  determined  value  would  be  highly  desirable. 

It  must  be  emphasized  that  the  instrument  error  model  used  assumes  linear  calibrations 
with  consunt  bias  and  scale  factor  errors  and  excludes  such  effects  as  non-linear  or  Mach  number 
dependent  calibrations,  friction  or  deadband  effeett  and  so  on.  While  these  assumptions  are 
probably  reasonable  for  the  manoeuvre  considered,  the  assumption  of  invariance  with  Mach 
number  of  airspeed  and  angle  of  attack  calibrations  may  be  doubtful  at  transonic  speeds, 
implying  that  speed  ranges  should  be  kept  to  a  minimum,  a  condition  unfavourable  to  lag 
determination.  Thus  regions  of  non-linear  calibrations  should  be  avoided  as  far  as  possible  to 
ensure  accuracy  of  the  procedure. 

The  question  of  the  accuracy  of  the  extracted  parameters  has  only  been  treated  briefly  and 
this  important  question,  including  the  effect  of  lags  on  the  parameter  values,  will  be  addressed 
more  fully  in  a  subsequent  report. 
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APPENDIX  1 
ScwMvity  Matrix  (SENS) 


The  parameter  vector  has  dimension  18  (equation  (S)): 

?  =  [A«x,  ft*x»  h»i»  bn.  Ay,  bv.  A,,  bg,  A^  bg,  Ain  iKo),  i^o),  8(o),  Xa,  Za]^ 
and  the  output  vector  has  dimension  4  (equation  (4)): 

i  =  [1^0Ut>  >OUt<  ^OUt»  AoutP 

It  follows  that  the  4 X  18  sensitivity  matrix  A  is  defined  by: 
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and 


/l  =  {»V-K1  +Aq)fll„  +  6,lx,}/M 


1^0Ut.{j  =  (1  +  AvX«  .U(j  +  W.  + 


'«o««.{j  =  (1  +^.XH>{j-/l  .  M{j)Ml  +/l*) 
0oat.{j  =  (1  +  A#)tf{j 


Aont.fj  =  (I  +A||)h|j. 


These  relationships  hold  in  general  for  those  elements  (j  of  the  parameter  vector  as  listed 
above  in  the  sensitivity  matrix.  The  only  exceptions  are: 


=  (I  +  AaX**’*,-/!  •  «A,-4in  •  X,)lu{l  +/i*) 

aoot.Dq  =  (1  +  A,XH-ba-/l  .  «b,-x.)/«(l  +/l*). 


APPENDIX  2 

DMKirMtW  Ef— rial  (DERIVS) 

In  order  to  calculate  the  state  (equation  (I))  and  the  elements  of  the  sensitivity  matrix 
(Appendix  I)  the  following  system  of  39  simultaneous  ordinary  differential  equations  needs  to 
be  integrated  with  respect  to  time; 
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fz  ~  (1  +'^<l)0ni +Aa- 

Initial  conditions  are  zero  for  all  variables  with  the  exception  of  the  states  u(o),  h>(o),  0(o) 
which  are  among  the  parameters  to  be  estimated,  and  the  two  variables  uoio)  and  wit(0)  which 
are  initially  I  0. 
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16.  Abstract 

-This  note  considers  matters  related  to  the  application  of  instrument  compatibility  checking 
techniques  to  flight  test  data.  A  previously  developed  Maximum  Likelihood  program  has  been 
used  to  study  the  effects  of  the  presence  of  scale  errors,  accelerometer  offsets  and  measurement 
time  lags  using  simulated  data.  Some  additional  information  on  the  effects  of  noise  levels  has 
also  been  obtained.  The  results  have  led  to  a  suggested  method  for  determination  of  centre  of 
gravity  location  from  flight  data.  The  effects  of  measurement  lags  have  been  shown  to  have  a 
major  influence  on  extracted  instrument  parameters  and  a  systematic  procedure  for  the  deter¬ 
mination  of  relative  phases  has  been  devised  and  applied  successfully  to  simuhted  data.  These 
techniques  have  also  been  applied  to  flight  data  from  a  roller-coaster  manoeuvre  and  a  set  of  rela¬ 
tive  lag  values  clearly  identified.  The  question  of  the  accuracies  of  extracted  instrument  parameters 
and  their  dependence  on  the  relative  lags  will  be  treated  more  fully  in  a  stdtsequent  publication. 
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